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ABSTRACT 

The surface current response to winds is analyzed in a two-year time series of a 12 MHz (HF) Wellen 
Radar (WERA) off the West coast of France. Consistent with previous observations, the measured 
currents, after filtering tides, are of the order of 1.0 to 1.8% of the wind speed, in a direction 10 to 40 
degrees to the right of the wind, with systematic trends as a function of wind speed. This Lagrangian 
current can be decomposed as the vector sum of a quasi-Eulerian current Ue, representative of the 
top 1 m of the water column, and part of the wave-induced Stokes drift U ss at the sea surface. Here 
Uss is estimated with an accurate numerical wave model, thanks to a novel parameterization of wave 
dissipation processes. Using both observed and modelled wave spectra, U sa is found to be very well 
approximated by a simple function of the wind speed and significant wave height, generally increasing 
quadratically with the wind speed. Focusing on a site located 100 km from the mainland, the wave 
induced contribution of U S s to the radar measurement has an estimated magnitude of 0.6 to 1.3% of 
the wind speed, in the wind direction, a fraction that increases with wind speed. The difference Ue of 
Lagrangian and Stokes contributions is found to be of the order of 0.4 to 0.8% of the wind speed, and 45 
to 70 degrees to the right of the wind. This relatively weak quasi-Eulerian current with a large deflection 
angle is interpreted as evidence of strong near-surface mixing, likely related to breaking waves and/or 
Langmuir circulations. Summer stratification tends to increase the Ue response by up to a factor 2, 
and further increases the deflection angle of Ue by 5 to 10 degrees. At locations closer to coast, 
Uss is smaller, and Ue is larger with a smaller deflection angle. These results would be transposable 
to the world ocean if the relative part of geostrophic currents in Ue were weak, which is expected. 
This decomposition into Stokes drift and quasi-Eulerian current is most important for the estimation of 
energy fluxes to the Ekman layer. 



1. Introduction 

D Surface drift constitutes one of the most important ap- 
plications o f the emerging oper ational oceanography sys- 
tems (e.g. lHackett et al. 2006), as it plays an important 
role in the fate of oil pollutions and larvae recruitment. A 
quantitative understanding of the relative contribution of 
the wave-induced Stokes drift to the near surface veloci- 
ties is also para mount for the prop er estimation of air-sea 
energy fluxes (Kantha et al. 2009). The quantitative vari- 
ation of surface drift as a function of the forcing parame- 
ters is still relatively poorly known. In areas of strong cur- 
rents due to tides or quasi-geostrophic dynamics, the sur- 
face drift current is highly correlated to the sub-surface 
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current. Otherwise, winds play a major role in defining 
the surface velocities. 

Re c ent theoretical and numerical works (|Ardhuin et"ai] 
[2001: iKantha and Clavsonl 12001: iRascle et all 120061: 
lArdhuin et alJ l2008bl) have sought to reconcile histor- 
ical measurements of Eulerian and Lagrangian (i.e. 
drift) vel ocities with recent kn owledge on wave-induced 
mixing ([Agrawal et al.l [1992) and wave-induced drift 
dRascle et al.l 120081) . These suggest that the surface 
Stokes drift U ss induced by waves typically accounts 
for 2/3 of the surface wind-induced drift, in the open 
ocean, and that the surface wind-related Lagrangian ve- 
locity Ul{z) is the sum of the strongly sheared Stokes 
drift Us{z) and a rel atively uniform quasi-Eulerian cur- 
rent u(z), defined bv l JenkinsI (1 19871) and generalized by 
lArdhuin et al.l d2008bl) . The Stokes drift decays rapidly 
away from the surface on a scale which is the Stokes 
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depth D$. For deep-water monochromatic waves of 
wavelength L, we take D$ = L/4, by analogy with the 
usual definition of the (twice la rger) depth of w ave influ- 
ence for the orbital motion (e. g. lKinsmanl 1 9651) . Namely, 
at that depth, the Stokes drift is reduced to 4% of its 
surface value. For random waves, a similar result re- 
quires a more complex definition, but the approximate 
same result can be obtained by using the mean wave- 
length L03 = gT? n03 where T m o3 is the mean period de- 
fined from the third moment of the wave frequency spec- 
trum (s ee Appendix). Sma ller values, like L/(4tt) used 
by e.g. iPolton et al.l d2005l) . are more reprensentative of 
the depth where the Stokes drift is truly significant. 

For horizontally homogeneous conditions, the depth- 
integrated quasi-Eulerian mass transport vector M m is 
constrained by the balance between th e Coriolis force and 
the wind (r a ) and bottom (77,) stre sses dHasselmannl 1 970t 
lArdhuin et al1l2004l ISrnkhll2006h . 



dM r > 



+ (M m + M w ) xe z =r a 



(1) 



where M™ is the (Stokes) mass 'transport^ induced by 
surface gravity waves, / is twice the vertical component 
of the Earth rotation vector, usually called the 'Coriolis 
parameter', and e z is the vertical unit vector, pointing up. 
The surface stress vector r a is typically of the order of 
p a CdU? with p a the air density and Cd in the range 1- 
2x 10 and Uw the wind speed at 10 m height. The hor- 
izontal homogenei ty is obviously never achieved strictly 
(e.g. iPollardl 1 1 9 8 3b . and this aspect will be further dis- 
cussed in the context of our measurements. 

The wind-driven current is not expected to be signif- 
icant at a depth greater than 0.7 times the Ekman depth 
D E = 0Aif(T a /p w )/ f (i.e. less than 0. 2% of the win d 
speed if the surface value is 2.8% of £/ in . lMadsen|[T977l) . 
For a wind speed Uio = 10 m s _1 , 0.7De is of the order 
of 30 m. In locations with a larger water depth, the bot- 
tom stress is thus expected to be negligible. Further, this 
depth of maximum influence can also be limited by a ver- 
tical stratification, with larger velocities in shallow mixed 
layers, and directions of Ue more strongly deflected 
to the right of the wind (in the Northern Hemisphere) 
than previou sly expected (Price and Sundermeve rll 1 999t 
lRascldl2007l) . It has also been proposed bv lPolton et al.l 
(2005) that the wave-induced mass 'transport' M w may 
play a role in the modification of near-surface currents, 
but M w is generally less than 30% of the Ekman trans- 
port M E = T a j f, and its effe ct appears to be secondar y 
compared to the stratification (Rasc le and Ardhuinl2009l) . 
The time-averaged balance given by (1) is thus approxi- 
mately, M m = -M w + (r a x e z ) If. This was nearly 
verifi ed for the LOTUS 3 dataset (Pric e and Sundermeved 
1999), when allowing for wave-induced biases in the 



'Because in the momentum balance (T} the term M™ drives a com- 
ponent of mean transport that opposes M™, there is no net wave- 
induced transport, except i n non-stationary or no n-homogenous con- 
ditions iHassehnanrJ )1970l) ; IXu and BoweiiUl994l) . 



mooring measurements dRascle and Ardhuinll2009l). Yet 
this i s not always the case (e.g. lNerheim and Stigebrandtl 
2006), possibly due to baroclinic currents and other phe- 
nomena that are difficult to separate from the wind-driven 
component. 

The vertical profile of the quasi-Eulerian current is, un- 
der the same homog eneous and sta t ionary circumstances, 
the so lution of (iXu and Bowed 11994 lArdhuin et al.l 
2008B 



du 

~dt 



+ (u + u s ) x e z 



d_ 

dz 



K 



dz 



(2) 



where K is a turbulent mixing coeffi cient. 

These predictions were verified bv lRasclel d2007) with 
mooring data at depths gr eater than 5 m and surfa ce- 
following measurements by Santala and Terravl (fl992) at 
depths larger than 2 m. When extrapolated to the sur- 
face using a simple numerical model, these observations 
give directions of Ue between 45° and 90°, more t han the 
45° gi ven by the const ant eddy viscosi ty model of[Ekman 
d 19051) . as extended by iGonellal d 1 97 lb . and the 10° give n 
by the linear eddy viscosity model of iMadsenl d!977l) . 
This surface angle, and the magnitude of Ue is also criti- 
cal for the estim ation of the flux of win d energy to the Ek- 
man layer (e.g. lWang and Huang 2004), or the analysis of 
near-surface drifter data ( e.g. iRio and Hernandez! 120031: 
Elip ot and Lumpkinl2008l) . For a better understanding of 
these questions, it is thus necessary to use ocean veloci- 
ties measured much closer to the surface. 

High Frequency (HF) radars can provide such mea- 
surements, at depths that depen d on their operating fre- 
quency. Using a 30 MHz radar. iMao and He ron (2008) 
made observations that are also consistent with the idea 
that the drift current, found to be 2. 1 % of the wind speed 
on average, is the sum of Ue which, according to their 
theory, depends quadratic ally on the wind speed, and 
U ss which they estimate to depend linearly on the wind 
speed, with a variation according to the fetch. Unfortu- 
nately, their analysis relied on empirical wave estimates 
that give large relative errors (of the order of 100% , see 
e.g. Ikahma and Calkoenl 199a lArdhuin et alJ2007h . and 
a limited range of wind speeds. Other HF-radar observa- 
tions gi ve a surface current of the order of 1.5 to 2.5% 
of U w dEssenll 1993b with 25 to 30 MHz radars. Dobson 
et al. (19 89) also report a ra tio of 2.0% using a 22 MHz 
radar, and lShav et al.l d2007l) report a ratio of 2 to 3% us- 
ing a 16 MHz radar in water depths of 20 to 50 m. These 
analyses are difficult to interpret due to the filters applied 
on time series to remove motions (tides, geostrophic cur- 
rents ...) that are not related to the wind, and also because 
of the importance of inertial oscillations that make the 
wind- and wave-driven current a function of the full wind 
history, and not just a function of the wind vector at the 
same time and location. 

In the present paper we extend the previous analyses 
of HF radar data by independently estimating the Stokes 
drift, using an accurate wave model. We find that at our 
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deep watei0 North-East Atlantic site the quasi-Eulerian 
current U E is of the order of 0.6% of the wind speed with 
a direction that is, on average, 60° to the right of the wind. 
We also find that the time-dependent response of surface 
current to the wind is typical of a slab layer with a trans- 
fer function proportional to l/(/ + where u> is the 
radian frequency considered. This result is expected to 
be representative of the open ocean. Therefore the es- 
timates o f the flux of wind energ y to the Ekman layer 
by e. g. IWang and Huang! (|2004) may not be quantita- 
tively correct: they used an angle of 45°, a surface veloc- 
ity which is 2y /, r a /p w for steady winds (about 0.2% of 
the wind speed), and a transfer function proportional to 
l/y/f + ijJ. A proper analysis of the effects of waves is 
needed to properly evaluate energy fluxes. 

Our new data and its processing are described in sec- 
tion 2, and the analysis of the stratification effect is pre- 
sented in section 3 with conclusions in section 4. 

2. Lagrangian and quasi-Eulerian current from HF 
radars 

a. Radar measurements and processing 

Hi gh frequency rada rs measure, among other things 
(e.g. Ilvonin et al.ll2.004y) . the phase velocity C of Bragg 
waves that have a wavelength equal to one half of the 
radar electromagnetic wavelength and that propagate in 
directions away from and toward the radar. This phase ve- 
locity is a combination of the quasi-Eulerian current Ue 
(IStewart and Joyll 1 974t iKirby and Chenlll989t) . the phase 
spee d of linear waves Cij n , an d a nonlinear wave correc- 
tion dWeber and Barrickll 19771) that can be interpreted as 
a filtered surface Stokes drift Usf- For monostatic sys- 
tems, the usual radial current velocity in the direction 9b 
towards one radar can be expressed as 

U R {9 B ) = C(6 B )-C rin -e 0B 

= U sf (e B ) + U E -ee B , (3) 

where eg B is the unit vector in direction 6b- This velocity 
can be loosely interpreted as the projection in direction 
9b of a current vector Ur. The reason why this is not 
exactly true is that Us /(0b) for all directions cannot be 
exactly given by the projection of a vector Usf- In other 
words, Us /(Ob) is not exactly proportiona l to cqs(0b), 
although it is a reasonable approximation dBroche et al.l 
119831 

In order to express Usf, we first define the Stokes drift 
vector for waves with frequencies up to f c from the di- 
rectional wave spectrum E(f, 8), 

U ss (f c ) = in [' P fk(f, 9)E{f, 9)dfdd, (4) 
Jo Jo 

where k(f) is the magnitude of the wavenumber k, equal 
to (2wf) 2 /g for linear waves in deep water, and g is the 

2 This means deeper than both the Stokes depth D 5 and the expected 
Ekman depth D e ■ 
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accelerat ion of gravity. Starting from the full exp r ession 
given bv I Weber and Barrick! d 19771) . iBroche et all d 1983b 
showed that the filtered Stokes drift component that af- 
fects the radial current measured by one radar station is 
well approximated by 

Usf(kB,0 B ) - U ss (f B )-e 0B 

poo p2ir 

+ 47rfc B / / fcos(e-6 B )E{f,e)d6df 
Jj B Jo 

(5) 

where fg is the frequency of the Bragg waves, and ks 
is the corresponding wavenumber vector, with a direction 
9b and magnitude k s- The full expres sion, correcting 
typographic errors in lBroche et al.l d 19831) is given in Ap- 
pendix A. In order to simplify the notations, the variable 
kt in Usf will now be omitted, but the filtered Stokes 
drift is always a function of the Bragg wavenumber, thus 
being different for different radar frequencies. 

The depth-varying quasi-Eulerian current u(z) is de- 
fined as the difference of t he Lagrangian velocity and 
Stokes drift (Jenkins! 1 19871) . and can generally be es- 
timated from the full v elocity field using a G eneral- 
ized Lagrangian Mean (Ard huin et al.l l2008bl) . The 
value U E estimated from the radar is, according to lin- 
ear wave theory, the integral of u(z) weighted by the 
Bragg wave Stokes d rift profile dStewart and JovHl974t 
iKirbv and Chen|[T989l) . In deep water this is, 

,o 

U E = 2k B eg B ■ / ue 2ksZ dz. (6) 

J — oo 

Here we use da ta from a WERA HF-radar system 
dGurgel et al.ll 19991) . manufactured by Helzel GmbH, and 
operated at 12.4 MHz. The Bragg wavelength is 12.1 m, 
corresponding to a wave frequency of 0.36 Hz in deep 
water. Thus half of the weight c 2fcfsZ in eq. (0 comes 
from water depths less than 0.6 m from the moving sea 
surface, compared to .28 m with the 30 MHz radar of 
iMao and Heronl d2008l) . The relative contributions from 
deeper layers to Ue decrease exponentially with depth 
as exp(2fcsz). Therefore Ue can be interpreted as the 
quasi-Eulerian current in the top 1 m of the ocean. 

The radar system has been deployed and operated 
by Actimar SAS, since July 2006 on the west coast of 
France (figure 1), measuring surface currents and sea 
states every 20 minutes. The area is characterized by 
intense tidal currents, in particular between the largest 
islands where it exceeds 3 m s _1 during mean spring 
tides. Also important, the offshore stratification is largely 
suppressed by mixing due to the currents in the areas 
shallower than 90 m, resulting in complex temperature 
fronts that are related to t he bottom topography (e.g. 
Mariette and Le CanSCHl. 

Each radar station transmits a chirped continuous wave 
with a repetition frequency of 4 Hz and a 100 kHz band- 
width which gives a radial resolution of 1.5 km. The re- 
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ceiving antennas are 16-element linear arrays with a spac- 
ing of 10 m, giving a typical angular resolution of 15 de- 
grees. The raw dat a is processed to remove m ost of the in- 
terference signals (Gurgel and Barbin 2008). Ensemble- 
averaging over 4 consecutive segments of 512 pulses 
yields a velocity resolution d u = 0.09 m/s in the Doppler 
spectrum used to estimate each individual radial cur- 
rent measurement. Yet, the current value is obtained by 
a weighted sum over a 9-point window applied to the 
Doppler spectrum. Provided that some inhomogeneity 
exists in the current field, the width of the Doppler spec- 
trum permits a measurement resolution that is infinitely 
small, but with an accuracy that is difficult to define, be- 
cause no ot her instrumen t, except maybe for the CODE- 
type drifter dDavis|[l9 85), is able to measure surface cur- 
rent in the top one meter of the ocean. Similarly, satel- 
lite altimeters are reported to measure the mean sea level 
position with an accuracy of the order of 2 cm whereas 
their t ypical range resolution is close to 40 cm. iPrandld 
(1987) used the coherence of the tidal motions to infer 
that the accuracy of his 27 MHz radar system was indeed 
less than the Doppler resolution when averaged over one 
hour. We will thus take the accuracy to be equal to the 
resolution, but as it will appear below, the only source of 
concern for our analysis is not so much the random error 
but a systematic bias, since we will average a very large 
number of independent measurements. 

Because we investigate the relationship between sur- 
face currents and winds based on modelled winds and 
waves, we will consider only the temporal evolution of 
the wave field at one point of the radars' field of view 
that is representative of the offshore conditions, at a dis- 
tance of 80 to 100 km from shore and with a water depth 
of 120 m. The reason for chosing this location is that we 
have verifed the wind and wave model results to be most 
accurate offshore where they were verified in situ with 
measurements that only span 6 and 9 months of our radar 
time series. Other reasons for looking at offshore condi- 
tions are the expected limited effect of the bottom, and 
the expected small horizontal gradients of both tidal cur- 
rents and other processes. Namely, we stay away from the 
therm al front that typically follows the 90 m depth con - 
tour dMariette and Le CanrJ[i985t iLe Bover et~al]|2009l) . 
The down side of this choice is that the HF-derived cur- 
rent is generally less accurate as the distance from the 
coast increases, and the coverage is not permanent, es- 
pecially during severe storms (e.g. figure Q]). These two 
drawbacks are limited in practice, as we now discuss. 

Interferences and ships cause some data to be rejected 
in the radar processing, or yield bad measurements, and 
heavy seas or calm seas also reduce the working radar 
range. In order to obtain a nearly continuous time series, 
we compiled and filtered data from a 0.2° in latitude by 
0.3° in longitude box around that point (A in figure 1, the 
arrow spacing indicate the resolution of the radar grid). 
This compilation was done in two steps. First, based on 
a visual inspection of the data, at each radar grid point, 
0.05% of the total number of data points in the radial ve- 



locities time-series are considered spurious and removed. 
These points are selected as the points where the raw ra- 
dial current time-series differs most from the result of a 
5-point median filter. The 0.05% value was selected as 
a convenient rule-of-thumb, which removes most of the 
visibly spurious points, but does not introduce too many 
unnecessary gaps in the time-series. Second, the time- 
series of all the grid points in the box around A were con- 
verted to u and v components and averaged. 

The Cartesian components of Ur and U e with respect 
to west-east (component u) and south-north (v) directions 
are calculated from the two radial components Ur(0bi) 
and Ur(0b2), each measured by one radar station, be- 
fore and after the substraction of U sf (&b ) ■ These Carte- 
sian components suffer from a geometr ical dilution of 
precision (GDOP), vary ing with position ( Chapm an et al.l 
119971 : IShav et aLTl2007l) . The radar beams intersect at 
point A with an angle r = 34° and it is possible to 
estimate the GDOP values for u and v, i.e. the ra- 
tios S u /s and S v /s where S u , S v and S are the uncer- 
tainties in u, v and u r , respectively. Assuming that S 
has no bias and is uniformly distributed from —d u /2 to 
+d u /2, each radar measurement has an intrinsic uncer- 
tainty S u = 0.04 m s _1 and S v = 0.11 m s _1 . 

This compiled time series, extending from July 5 2006 
to July 31 2008, is the basis of the following analysis. 
The 1200 s resolution data was averaged over 3 h blocks 
centered on round hours. Gaps shorter than 6 h were lin- 
early interpolated. That time series is 97% complete, and 
thus covers two full years. Other parts of the radar field of 
view yield similar results, briefly discussed below. Due 
to the averaging in space and time, each point in the time 
series is the combination of about 30 range cells and 9 
time intervals, i.e. 180 independent velocity measure- 
ments when the full radar range is obtained. Even with 
all cm s _1 uncertainty on the original measurement, the 
expected r.m.s. error on the velocity components are thus 
less than 1 cm s _1 . This analysis assumes that the in- 
strument is not biased. After verification of the radar an- 
tenna lobe patterns using both in situ transmitters and a 
novel technique based on the analysis of radio interfer- 
ence (to be described elsewhere), the main lobe of the 
radar is known to be mispointed by less than 5 degrees, 
with a -3dB width less than 15°. The largest source of un- 
certainty is thus the interpretation of the phase speed and 
the numerical estimation of the Stokes drift, as discussed 
below. 

Because we wish to focus on the random wind-driven 
currents, we als o performed a tidal ana lysis using the T- 
TIDE software (Pa wlowicz et a l. 2002) applied to each 
velocity component. This analysis on the full time series 
(before time averaging) allows the removal of the deter- 
ministic diurnal constituents K\, Ox, Pi and Qi that have 
amplitudes of 1.5 to 0.3 cm s _1 , with estimated errors of 
0.1 cm s _1 . Because this only corrects for 95% of the 
apparent variance in the M% and S2 semi-diurnal tides, 
these will be further filtered using a time filter. 
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FIG. 1. Map of the area showing a map of significant wave height 
on January 1st 2008, at 12:00 UTC, estimated with a numerical wave 
model (see Appendix B), and the instantaneous surface current mea- 
sured by the H.F. radars installed at Porspoder and Cleden-Cap-Sizun. 
In situ measurement stations include the weather buoy BEAtrice and the 
Pierre Noires (62069) directional Datawell waverider buoy (installed 
from November 2005 to March 2006 and back again since January 
2008), and a previous waverider deployment (Iroise), more represen- 
tative of the offshore wave conditions. The large black square around 
point A is the area over which the radar data has been compiled to pro- 
vide the time series analyzed here, representative of offshore conditions. 
When the radar functionned, over the entire square measurements are 
available for more than 80% of the 20 minute records, a number than 
rises to 99% for the area East of 5°35'W. The partial radar coverage 
around point A is typical of high sea states with H s > 6 m offshore, 
which are rare events. 



b. Numerical wave model and estimations of Stokes drift 



1) General principles 

As expressed by eq. (01, the estimation of Us /(Ob) 
requires the measurement or modelling of the wave spec- 
trum E(f, 6). In situ buoys were moored for restricted 
periods at several locations for the i nvestigation of off- 
shore to coastal wave transformation ( Ardhuin 2003) and 
to provide complementary data for radar validation. The 
radar also measures the sea state, but the coverage is often 
limited, and its accuracy for a 20 minute record is typi- 
cally only of the order of 25% for the significant wave 
height H s . Thus, in order to use the full current time se- 



ries at the offshore location (point A) we have to estimate 
the sea state using a numerical wave model. 

We use an implemen tation of the WAVE WATCH III 
code, in its version 3.14 jTolmanl2007l l2008). with minor 
modifications of the parameterizations, see appendix B, 
and the addition of advection schemes on unstructured 
grids (lRolandll2008l) . 

The model setting consists of a two-way nested pair of 
grids, covering the global ocean at 0.5 degree resolution 
and the Bay of Biscay and English channel at a resolu- 
tion of 0.1 degree. A further zoom over the measurement 
area is done using an unstructured grid with 8429 wet 
points (figure 1). The model setting is fully described in 
appendix B. 

In practice, Usf is dominated by the first term 
Uss(fB), in eq. (0. Examining a large number of spec- 
tral data (6 buoys for 2 years spanning a range of wave 
climates, see appendix C), we realized that C/ ss (/s) is es- 
sentially a function of the wind speed U±o and the wave 
height H s . While U\q explains typically only 50% of the 
variance of U ss (f) with 0.3 < / < 0.5, U±o and H s gen- 
erally explain over 85% of the variance. This behaviour 
of U ss (f) is similar to that of the fourth sp ectral moment, 
related to the surface mean s quare slope dGourrion et alj 
2002; Vandem arket1d1l2004]) . The reason for this corre- 
lation is that the wind speed is obviously related to the 
high frequency part of the wave spectrum, which deter- 
mines most of the Stokes drift, while H s is a surrogate 
variable for both the presence of swell and the stage of 
development of the wind sea. Here we find, 



U ss (f c ) 



5.0 x 10" 4 



1 .2.'. 0.25 | ^ 

Jc 



1.3 



u 



10 



min {U 10 , 14.5} + 0.025 (H s - 0.4) . 



(7) 

The relationship given by eq. (Q appears to be very ro- 
bust, with a 2.6 cm _1 r. m. s. difference compared to 
global hindcast values of U ss (oo), which is a 16.9% dif- 
ference. Nevertheless, when compared to buoy data, an 
accurate wave model generally provides a better fit to the 
observations (Appendix C). We thus have used our hind- 
casts using WAVEWATCH III to provide an estimate for 
Usf- 

2) Uncertainty on Usf around point A 

We have no wave measurement at point A, and no per- 
manent spectral measurement in the area. A detailed val- 
idation of U ss was thus performed for the coastal buoys 
62069 (figure 1), 62064 (off Cap Ferret, 600 km to the 
southeast of point A), the U.S. Northwest Pacific Coast 
(appendix C), U.S. East coast, Gulf of Mexico and Cali- 
fornia. 

We further use wave information at buoy 62163, lo- 
cated 150 km west of point A, representative of the off- 
shore conditions found at point A, and a combination of 
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satellite altimeter data. The present model estimates of 
H s are more accurate at buoy 62163, located 150 km 
west of point A, than at Pacific buoy locations. Fur- 
ther, the model estimate of the fourth moment m 4 of the 
wave spectrum is better correlated in the Bay of Biscay to 
radar altimeter C-band cross-section, compared to other 
regions of the world ocean (Appendix C). We thus expect 
the model estimate of U ss (fB = 0.36 Hz) to have a bias 
smaller than than 5%, with a random error less than 20% 
(see Appendix C). As a result, We chose to use this nu- 
merical wave model for the estimation of U ss and Usf. 
We can thus propose an error buget for our estimate of 
the wind-driven quasi-Eulerian current in which the mea- 
surement error is dominated by Usf with a bias of 5% 
at most and a standard deviation less than 20% overall. 
Using the analysis of 2 years of model results, this stan- 
dard deviation at the Pacific buoy 46005 is 24% for wind 
speeds of 3 m s -1 , 20% for 5 m s _1 , 16% for 7 m s~\ 
11% for 11 m s _1 . Given the general accuracy of the 
wave model in the North-East Atlantic, we expect similar 
results here. 

We thus estimate that the root mean square error of the 
modelled quasi-Eulerian current Ub at 3 hour intervals 
is of the order of 0.2% of U\q. On this time scale, it is 
difficult to rule out contributions from horizontal pressure 
gradients in the momentum balance, and this current may 
not be purely wind-driven. 

The averaged current, e.g. for a given class of wind 
speed, as shown on figure 7, has a relative accuracy bet- 
ter than 0.1% of U\q. In-situ measurements of time- 
averaged velocities from 10 to 70 m above the bottom 
at 48°6'N and 5°23'W (south of point A, see figure 1) 
using a RDI Wo rkhorse ADCP deplo yed from June to 
September 2007 dLe Boyer et al.ll2009l) give tide-filtered 
currents less than 2 cm s _1 or 0.25% of the wind speed 
when averaged following the wind direction (the instan- 
taneous measurements are rotated before averaging), and 
less than 0.1% when winds stronger than 10 m s -1 . 
This is typically less than 20% of Usf- Assuming that 
wind-correlated baroclinic currents are negligible during 
the ADCP measurement campaign, the wind-correlated 
geostrophic current is expected to be less than 0.2% of 
Uio. Gereralizing this result to the entire radar time se- 
ries, the averaged values of Ue can be interpreted as a 
wind-driven current with an accuracy to within 0.3% of 
Uio. 

3. Analysis of wind-driven flows 

The study area is characterized dominated by moder- 
ate 6 to 12 m s _1 winds, from a wide range of direc- 
tions, with slightly dominant South- Westerly and North- 
Easterly sectors (figure |2). 

a. Rotary spectral analysis 

The rotary spectral analysis gives both the frequency 
distribution of th e signal, and an indication of its circu- 
lar polarization dGonellalll971l) . The positive frequen- 




FIG. 2. Wind rose for the years 2006 to 2008 at point A, based on 
ECMWF analyses. The observations at BEAtrice buoy give a similar 
result. For each direction, the cumulative frequency is indicated with 
wind speeds increasing from the center to the outside, with a maximum 
of 4.3% maximum from West-South- West (heading 250°). An isotropic 
distribution would have a maximum of 2.7%. 



cies correspond to counter-clockwise motions, and the 
negative frequencies correspond to clockwise motions, 
the usual polarization of inertial motions in the Northern 
Hemisphere. 

The instantaneous measurements of the radar are dom- 
inated by tidal currents, and the variance of motions with 
frequencies less than 1.75 count per day (cpd) only ac- 
counts for 8% of the total variance (figure 3). These low 
frequency motions include the diurnal tidal constituents, 
most importantly K\ and Ox, but these only account for 
0.1% of the variance. The low frequency motions are 
generally dominated by near-inertial motions, which are 
polarized clockwise with frequencies close to the inertial 
frequency fj = 1.3 counts per day (c.p.d., see figure 3). 

b. Co-spectral analysis 

Here we investigate the relationship between measured 
currents, processed as described above, and winds, taken 
from 6-hourly wind analyses from ECMWF. These anal- 
yses were verified to give excellent correlation (r ~ 0.92) 
with the BEA buoy (WMO code 62052), which unfortu- 
nately malfunctionned during large periods of time. The 
wind and current data are thus completely independent. 
The wave model was forced by these same winds, and 
thus the high level of coherence between the predicted 
Stokes drift and the wind (figure 4) is not surprising. 

In order to isolate the wind-correlated dynamics from 
the shorter (tide) and longer (general circulation) time 
scales, we first perform a co-spectral analysis of the mea- 
sured currents with the wind, following the method of 
iGonellal dl971l) . In order to keep as much data as possi- 
ble between data gaps, the Fourier transforms are taken 
over 264 hours, which corresponds to 21 M 2 tidal cy- 
cles. The measured currents are significantly coherent 
with the wind vector over the range -1.75 to 1.75 cpd 
(figure 4). This coherence is generally reduced when the 
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Stokes component Us / is subtracted from the radar mea- 
surements. 




FIG. 3. Rotary power spectra of the current measured by the rada r, and 
the contribution U g * to the surface Stokes drift estimated via eq. iAU . 
Clockwise (CW) motions are shown with dashed lines and counter- 
clockwise motions are shown with solid lines. The spectra were esti- 
mated using half-overlapping segments 264 h long over the parts of the 
time series with no gaps. The number of degrees of freedom is taken 
to be the number of non-overlapping segments, i.e. 59, at the spectral 
resolution of 0.09 cpd, giving a relative error of 35% at the 95% con- 
fidence level. In the bottom panel the the tidal components have been 
filtered out, which clearly removes the diurnal peak However, the the 
semi-diurnal tides are only reduced by a factor 25, which is not enough 
compared to the magnitude of the near-intertial motions, and requires 
the use of an additional filter. This tide-filtered time series is used in all 
of the following. 



The radar-measured current vectors U r have stable di- 
rections relative to the wind, 20 to 40° to the right for 
/ > — //, given by their coherence phase (figure 4). The 
coherence phase of the Stokes drift increases with fre- 
quency. This pattern is typical of a time lag, that can 
be estimated to about 1 .5 hours, consistent with the rel- 
atively slow response of the wave field compared to the 
current. This is rather short compared to the time scale 
of wave development, but one should bear in mind that 
the Stokes drift is mostly due to short waves that respond 
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FIG. 4. Rotary co-spectra of the wind and wind stress with the radar- 
derived current, Stokes drift and Eulerian current, (a) magnitude and 
(b) phase. The number of degrees of freedom is 108 at the spectral 
resolution of 0.09 cpd. Coherence is significant at the 95% confidence 
level for a value of 0. 1 . Negative and positive frequencies are clockwise 
and counter-clockwise polarized motions, respectively. 



faster to the wind forcing than the dominant waves. Be- 
cause the wind preferentially turns clockwise, the Stokes 
drift is slightly to the left of the wind. The asymmetry 
in the phase of Us / for clockwise and counter-clockwise 
motions may be related to varying fetch when the wind 

turns. 

As expected from the theory by iGonellal dl972f) . the 
phase of the quasi-Eulerian current Ue jumps by about 
180° at the inertial frequency — //. In the frequency 
range from -1.2 to 0.2 cpd, that contains 40% of the non- 
tidal signal, Ue is at an angle between 45 and 60° to the 
right of the wind. This conclusion is not much altered 
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when one correlates the Eulerian current against the wind 
stress, which, for simplicity is estimated here with a con- 
stant drag coefficient, t = 1.3 x 10 _3 Uio?7io- One may 
argue that the theoretical filtering of the Stokes drift is not 
well validated. A lower bound on the estimate of Us / can 
be given by removing the contribution from waves shorter 
than the Bragg waves. This has very little impact on the 
estimation of Ue- 

The observed coherence phases of Ue and C/ 10 are sim- 
ilar to the values given by Gonella (1972, figure 6 ) , base d 
on the constant eddy-viscosity model of Ekman ( 1905), 
but for the current considered at a depth as large as 25% 
of the Ekman depth. Since the radar measurements are 
representative of the upper 1 meter, and the Ekman depth 
is generally of the order of 30 m, it follows that the 
classical Ekman theory, with a constant eddy viscosity, 
does not apply here. Instead, this large near-surface 
deflection is consistent with model results obtained with 
a high surf ace mixing such as induced by Langmuir cir- 



culati ons dMcWilliams et alj Il997h iKantha and Clavsonl 
20041). breaking waves ([Craig and Banner 1994; 



Mellor and Blumberg 2004; Rascle_etalJ|2006) or both, 
and consistent with the few observed near-surface 
velocity profiles (Santala and Terra vl[T992h . 

c. Effects of stratification 

Following the the ory of Gonella (1972) and the p revi- 
ous observations by Price and Sundermeyer ( 1999), it is 
expected that the stratification has a significant effect on 
the surface currents. Here we used sea surface tempera- 
ture time series to diagnose the presence of a stratifica- 
tion. Because of the strong vertical mixing year-round 
at the site of buoy 62069, the horizontal temperature dif- 
ference between points A and point 62069 is a good in- 
dicator of the vertical stratification at point A. This tem- 
perature difference reaches up to 2°C, and was present 
in 2006, 2007 and 2008 from early July to late October, 
as revealed by satellite SST data. We thus separated the 
data records used for the spectral analysis into "stratified" 
and "homogeneous" records based on the date of the mid- 
point in these time series. 

These two series show a significant difference (at the 
95% confidence level) when the spectra are smoothed 
over 0.3 c.p.d. bands, with a twice larger response in the 
cases expected to be stratified (dashed lines, figure 5) for 
frequencies in the range - 1 .7 to 1 .5 c.p.d. Interestingly the 
transfer functions decrease like l/(/ + ui) from a peak at 
the inertial frequency /, where u> is the radian frequency. 
This decrease is typical of slab-like behaviors that are ex- 
pected in mixed layers with a much larger surface mix- 
ing (e.g. Rascle et al. 2006) than typically used with 
Ekman theory, or a mixed layer depth much shallower 
than the Ekman depth (Gonella 1972). Ekman theory 
in unstratified conditions, that should apply to our win- 
ter and spring measurements, would give a much slower 
decrease, proportional to + lu) (Gonella 1972). 

Together with this stronger amplitude of the current re- 




FlG. 5. Amplitude transfer functions (top) and coherence phases (bot- 
tom) between the wind forcing and the current response. The dashed 
lines correspond to records where a stratification is expected to be im- 
portant (18 out of 108), and the solid lines correspond to the other 
records. Confidence intervals for the two group of records are shown 
for the native spectral resolution of 0.09 c.p.d. In order to be at a com- 
parable level the wind stress was multiplied by 50 before estimating the 
transfer function. The two peaks of the transfer functions at +/- 2 cpd 
are due to the tidal currents but do not correspond to a causal relation- 
ship between the wind forcing and the current response. 



sponse in stratified conditions, we find a larger deflection 
angle in the -0.8 to -0.2 c.p.d. frequency range. This 
pattern of larger currents and larger deflection angles in 
stratified conditions is consiste nt with the observations of 
iPrice ancLS undermever ( 1999), and t he numerical model 
results by Rascle and Ardhuir] d2009h . 
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d. Relationship between tide-filtered currents and winds 

A proper model for the wind-induced current may be 
given by the relationship between the wind speed and 
wave height, giving the Stokes drift, and the complex 
transfer function (transfer function and phase) from the 
wind stress spectrum to the Eulerian current spectrum , 
following Gonella (1971) or iMillot and Creponl (1198 ll) . 
Such a model is beyond the scope of the present paper. 

Simpler models that would give the current speed and 
direction as a function of the instantaneous wind vector 
are even less accurate. Because the transfer function is 
very peaked at the inertial frequency, the current speed 
may vary widely for a given wind speed. Yet, for prac- 
tical reasons, there is a long tradition of directly com- 
paring current and wind magnitudes and directions for 
search and rescue operations and ocean engineering ap- 
plications. Because of the inertial oscillations, there is 
usually a large scatter in the correlation of the current 
and wind speed vec tors. In order to com pare with pre- 
vious analyses (e.g. iMao and Heronlf2008l) . we thus per- 
form such a comparison, after filtering out the dominant 
tidal current, by taking the inverse Fourier transform of 
the current, wind, and Stokes drift spectra in which the 
amplitudes of components with frequencies higher than 
1.75 cpd, and the zero frequency, are set to zero. Again, 
the Fourier transforms are taken over 264 hours. 

We find that the surface Eulerian U e current lies 40 
to 60° to the right of the wind, suggesting that the near- 
inertial motions only add scatter to the longer period mo- 
tions (|/| < 1.3 c.p.d.) that were found to have sim- 
ilar deflection angles. Interestingly, the typical magni- 
tude of Ue decreases from about 0.8% of Uw at low 
wind to nearly 0.4% for high winds. This reduction in 
the relative magnitude of Ue is accompanied by a re- 
duction of the deflection angle from 65° on average for 
Uio = 3 m s" 1 to 40° for U 10 = 15 m s" 1 . On the 
contrary, the Stokes drift typically increases quadratically 
with the wind speed. These o bservations co ntradict the 
usual theoretica l state ments of iKirwan et al.l (1 19791) and 
IMao and Heronl d2008l) : they concluded that the Stokes 
drift should be linear and the Eulerian current should 
be quadratic in terms of wind speed. The fact that the 
Stokes drift is quadratic as a function of the wind speed 
is well shown by observed wave spect ra in figure CI 
and th e fitted equation (IT). The error in Ma o and Heronl 
(2008) is likely due to their erroneous assumption that 
the Stokes drift is dominated by waves at the pea k of 
the spectrum In th e analysis of IKirwan et al.l (1 19791) and 
iRascle et al.l (120061) . the error essentially arises from the 
assumed shape of the wave spectrum. 

The less-than-linear dependence of Ue on U\q con- 
tradicts the usual simple Ekman model for the quasi- 
Eulerian current, which would predict a current propor- 
tional to the wind stress, and thus varying as the square 
or cube of the wind speed. This difference is likely due 
to the enhanced mixing caused by breaking waves, which 
tends to mix the momentum over a scale of the order of 



the wi ndsea wave height, i.e. increasing wit h the wind 
speed dTerrav et all [19961 IRascle et al J 120061) . Numeri- 
cal models without stratification but with a realistic mix- 
ing tend to give a quasi-Eulerian current that increases 
with wind speed and with the inverse wave age. Here the 
stronger winds do not correspond to very different wave 
ages, and it is likely that a correlation of deeper mixed 
layers with stronger winds is the ca use of the reduction 
of Ue with increasing wind speed (IRascle and Ardhuinl 
120091) . As a result, the nonlinear current response to the 
wind stress will likely limit the accuracy of models based 
on transfer functions. 
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FIG. 6. Observed tide-filtered quasi-Eulerian velocity magnitudes, nor- 
malized by the wind speed, and directions, relative to the wind vec- 
tor. The linear increase of Us f /Uio with Uio is consistent with the 
quadratic dependence of Usf on Uio given by eq. (7J. The full dataset 
was binned according to wind speed. Dash-dotted lines correspond to 
stratified conditions only and dotted lines correspond to homogeneous 
conditions. The number of data records in each of these cases is indi- 
cated in the bottom panel. The dashed line show results when Us f is 
replaced by C/ss(/b). Error bars show only 1/2 of the standard devi- 
ation for all conditions combined, in order to make the plots readable. 
All time series (wind, current, Us f and U ss were filtered in the same 
manner for consistency (except for the initial de-tiding applied only to 
the current data). The error bars do not represent measurement errors 
but rather the geophysical variability due to inertial motions. 
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FIG. 7. Mean wind-correlated current vectors in low and high wind 
conditions, with and without stratification, measured off the West coast 
of France with the 12.4 MHz HF radar, based on the results shown in 
figure 6. Ur is the radar-measured vector, that can be interpreted as 
a sum of a quasi-Eulerian current U e > representative of the upper two 
meters, and a filtered surface Stokes drift U gt. The full surface Stokes 
drift is typically 40% larger that this filtered value. Solid circles give 
the expected error on the mean current components due to biases in the 
wave contribution to the radar measurement. The dashed circle show 
the expected error on the interpretation of Ue as a wind-driven current, 
based on the ADCP measurements at depth of 60 to 120 m, assuming 
that the baroclinic part of the geostrophic current is negligible. 



e. Effects of fetch or wave development 

The same analysis was also repeated for other points 
in the radar field of view. For example at point B (figure 
1), the radar data quality is generally better, but where the 
wave model may have a bias of about 10% on U ss , and 
the ECMWF wind field may be less accurate. Point B is 
relatively sheltered from Southerly, and North-westerly 
waves, and the fetch from the East is 40 km at most. If 
we assume that the winds are accurate at that site too, we 
find that the radar-derived current is weaker relative to the 
wind, with Ur/Uxo typically smaller by 0.2% point (i.e. 
a ~ 15% reduction) compared to point A. This appears to 
be due to a reduction in Us f, which is only partially com- 
pensated for by a small increase in Ue- This difference 
between A and B nearly vanishes when only Westerly 
wind situations are considered (defined by winds within 
60° from the Westerly direction). 

4. Conclusions 

Using a 2 year time series of HF radar data, and a novel 
numerical wave model that is shown to reproduce the ob- 
served variability of the surface Stokes drift with wind 
speed and wave height, we have analyzed the wind-driven 
surface current. When tidal currents are filtered out, the- 
ory predicts that the measured velocities are a superpo- 
sition of a filtered Stokes drift U$f and a quasi-Eulerian 



current U e- With our 12 MHz radar, Us f is estimated to 
be of the order of 0.5 to 1.3% of the wind speed, with a 
percentage that increases linearly with wind speed. These 
values are a function of the radar wavelengths and would 
be larger, by up to 20%, with higher frequency radars 
that give currents representative of a shallower surface 
layer. The other component Ue is found to be of the or- 
der of 0.6% of the wind speed, and lies, in our Northern 
Hemisphere, at an average 40 to 70 degrees to the right of 
the wind, with a large scatter due to inertial oscillations 
that may be we ll modelled using a Laplace transform of 
the wind stress (Bro che et al.ll9 83). This large deflection 
angle is robustly given by the coherence phase for clock- 
wise motions in the frequency range from to the inertial 
frequency. 

When instantaneous currents are compared to the 
wind, the magnitude of Ue appears to decrease with wind 
speed but it increases when a stronger stratification is ex- 
pected (figure 6). These surface observations correspond 
to currents in the depth range to 1.6 m, and confirm 
previous analysis of deeper subsurface mooring data. If 
wind-correlated geostrophic current are negligible in our 
measurements, the shape of the classical picture of the 
Ekman spiral is not correct, and the surface layer is much 
more slab-like than assumed in many analyses, proba- 
bly due to the large w ave-induced mixing at the surface 
(Agra wal et al.|[l992l) . These findings are summarized in 
figure 7. 

If we neglect the wind-correlated geostrophic currents, 
which we deem reasonable, and interpret Ue as being 
purely wind-driven, our observations of Ue/Uiq at point 
A are expected to be representative of the open ocean, 
whereas in coastal areas and small basins, a less devel- 
oped sea state will lead to a smaller value of Us f and 
a larger value of Ue, as we observe at point B. Such 
a generic relationship of Ue and Uio is very important 
for a proper estimation of the energy flux to the mixed 
layer. Besides, on top of the wind stress work on the 
Ekman current, this energy flux should be dominated by 
the dissipation of w ave energy induced by breaking (e.g. 
iRascle et al.ll2008h . Also, there is the depth-integrated 
Stokes-Coriolis force which is equal to the product of the 
depth-integrated Stokes transport M™ = p w J XJ s (z)dz, 
and the Coriolis parameter /. This force is smaller than 
the depth-integr ated Coriolis forc e by about a factor of 
3 on average (IRascle et al.l 120081) . but that may give a 
comparable work due to the smaller angle between that 
force and the quasi-Eulerian current u(z). The accu- 
rate estimation of the surface Stokes drift using a numer- 
ical wave model also opens the way for a more accu- 
rate interpretation of space-borne measurements of sur- 
face currents using Doppler methods, that are contami- 
nated by a Stokes-like com ponent amplified 10 times or 
more (Chapronetal. 2005|). 
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APPENDIX A Nonlinear correction for the wave 
dispersion relation in a random sea state 

Based on the l o west order approximate theory of 
I Weber and Barrickl (1 19771) for deep water waves with 
/ ~ 2TTy/gk, the nonlinear correction to the phase speed 
of components with wavenumber k B and direction 9 B , 
can be expressed as an integral over the wave spectrum . 
Defining x = k/k B and a = 6>-6> R . dBroche etalj 19831 
their eq. A2) give the following expression, 



u Sf (k B ,e B ) = 




F(x,a)E(f,9)d9df, 



where, correcting for typographic errors, and using y — 

x x l 2 = f / f B and a — cos a, 

F(x,a) = y {2a — y + 3xa} 

x {(ya-x) (a e + (l + ey) 2 )/2 

+ (1 + ey) (1 + exa + ey (x + ea) — a e )} , 

(Al) 

with 

a e = (l + x 2 + 2exa) 1/2 . (A2) 

These expressions give the correct figures in lBroche et al.l 
(119831) . For x < 1 one finds that F(x,0) = 4x 3/2 , and 
for x > 1, F(x,0) = \x x l 2 , as p re viously given by 
lLonguet-Higgins and Phillipsl (119621) iHuang andT ung 
(119761) and Barrick and Weber] (119771) . As commented 
bv iBroche et al.l d!983l) . F(x, a) ~ F(x,0)cosa, with 
the largest errors occurring for x — 1 where F(x, a) > 
F(x,0) cosa for \a\ < tt/3, which, in our case makes 
Usf larger by 2 to 5% than the approximation given by 
eq. ©. 

APPENDIX B Parameterization and numerical 
settings for the wave models 
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a. Parameterizations 

The implementation of the WAVEWATCH III model 
used here was ran with source functions Si„ , S„i and Sds 
parameterizing the wind input, nonlinear 4- wave interac- 
tions and whitecapping dissipation. An extra additional 
dissipation term Sdb is also included to enhance the dis- 
sipation due to wave break ing in shallow water, based on 
Batti es and Janssenl (119781) . 

The paramete ri zation for S n i is taken from 
Hasselman n et alj (Q985), with a minor reduction of 
the coupling coefficient from 2.78 x 10 7 to 2.5 x 10 7 . 
The parameterizatio ns for Sj n and Sds a re very similar 
the ones used by lArdhuin et al.1 d2008al) . with modifi- 
cations to furth er improve the h igh frequency part of 
the spectrum dFilipot et al.l 120081) . Namely, the white- 
capping dissipation is b ased on recent obse rvations of 
wave breaki ng statistics dBanner et al.l l2000b . and swell 
dissipation (Ar dhuin et al.l 12009b . These model settings 
give the best estimates so far of wave heights, peak 
and mean periods, but also of parameters related to the 
high frequency tail of the spectrum (appendix C). The 
present model result s are thus a signifi cant i mprovement 
over th e results of IB idiot etaf] (120071) and iRascle et al.1 
(2008). The physical and practical motivations for the 
parameterizations will be fully described elsewhere, and 
we only give here a description of their implementation. 
We only note for the interested users, that the parameter 
settings given here tend to produce larger negative 
bi asses on H s for H„ > 8 m than the parameterization 
by lBidlotetaLl d2007h . Better settings for H s in extreme 
waves would be s u = and C3 = 0.5 (see below), but 
this tends to give too large values of U ss , which is why 
we do not use these settings here. 

The parameterizati on of Sj n is take n from [Janssen 
dHHb as modified by Bidlot et al. ( 2007), with some fur- 
ther modifications for the high frequencies, and the addi- 
tion of a wind output term S ml t (or "negative wind input") 
based on the observations by Ard huin et alj d2009l) . The 
source term is thus 

Pw K V ^ / 

xcos 2 (0-8 u )<jF(f,9) + S out (f,0), 

(Bl) 

where /3 max is a (constant) non-dimensional growth pa- 
rameter, k is von Karman's constant, in the friction 
velocity in the air, C is the phase speed of the waves, a 
is the intrinsic frequency, equal to 27r/ in the absence of 
currents, and F (/, 9) is the frequency-directional spec- 
trum of the surface elevation variance. In the present im- 
plementation the air/water density ratio is constant. We 
define Z — log(^) where /i is given by Janssen (1991, 
eq. 16), corrected for intermediate water depths, so that 

Z = log(fei) + k/ [cos (9 - 9 U ) « + z a )] , (B2) 
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where z% is a roughness length modified by the wave- 
supported stress t w , and z a is a wave age tuning parame- 
ter. The effective roughness Z\ is implicitly defined by 



U 



10 
zo 

Zl 



/10 m 
— log 

K V Zl 



max <^ a — , 0.0020 
I 9 

y/l- t w /t 



(B3) 

(B4) 
(B5) 



where r is the wind stress magnitude, t w is the wave- 
supported fraction of the wind stress, Uio is the wind at 
10 m height and g is the acceleration of gravity. 

The maximum value of zq was added to reduce the un- 
realistic stresses at high winds that are otherwise given 
by the standard parameterization. This is equivalent to 
setting a maximum wind drag coefficient of 2.8 x 10 ~ 3 . 
This, together with the use of an effective friction veloc- 
ity instead of in (IB2b are the only changes to 
the general form of Janssen's (1991) wind input. That 
friction velocity is defined by 



K(/)) 2 = 



\ u l e e 



rf f 2n S in (f,9') 



o Jo 



(B6) 



Here the empirical factor s u = 1.0 adjusts the shelter- 
ing effect of short wave s by long waves adapted from 
IChen and Belcher! |2000), and helps to reduce the input at 
high frequency, without which a balance of source terms 
would not be possible (except with a very high dissipation 
as in Bidlot et al. 2007). This sheltering is also applied 
in the precomputed tables that gives the win d stress as a 
function of U w and t w /t dBidlot et al.ll2007l) . 

The wind output term , is identical to the one used 
by lArdhuin et al.1 (l2008al), bas ed on the satellite obser- 
vations of lArdhuin et al.l (120091) . with an adjustment to 
Pacific buoy data. Namely, defining the Reynolds num- 
ber Re= 4it or b<2 rb / Va, where w or b and a or b are the sig- 
nificant surface orbital velocity and displacement ampli- 
tudes, and v a is the air viscosity, we take, for Re< 10 5 



Sout (/) 9) 

and otherwise 

Sout (/, 0) - 
where 



-1.2^ hkV2^)F(f,6). (B7) 
o„, L J 



^{l6f e <T 2 u olh /g}F(f,0), (B8) 

Pw 



fe = 0.7 f e ,GM + [0.015 - 0.018 cos(0 - 9 U )} u*/u olbl 

(B9) 



where f e ,GM is the friction factor given by Grant and 
Madsen's (1979) theory for rough oscillatory boundary 
layers without a mean flow, using a roughness length ad- 
justed to 0.04 times the roughness for the wind. This 
gives a stronger dissipation for swells opposed to winds. 

The dissipation term is the sum of the saturation-based 
term of lArdhuin et al . (2008a) and a cumulative breaking 
term S ds c o fFilipotetalJ d2008). It thus takes the form 



S ds (f,6) =aC ds 0.25 



»*{tSP-m} 

{^-1,0} 



-0.75 

xF(f,6) + S dSiC (f 



(B10) 



where 



B'(f,9) 



6+80° 



k 3 cos 2 (6 - 9') F(f, 9')C g /{2ir)d9\ 



19-80° 

B(f) = max{B'(f,9),9€[0,2ir[}: 



(Bll) 
(B12) 



and B r = 0.0009 is a threshold for the onset of 
breaki ng cons istent with the observ ations of lBanner et al.l 
(2000) 



and 



Banner et al. (2002), as discussed by 
Babanin and van der Westhuvsen (2008), when including 
the normalization by the width of the directional spec- 
trum (here replaced by the cos 2 factor in eq. IB1 II ). 

The dissipation constant C ds was adjusted to 2.2 x 
10~ 4 in order to reproduce the directi onal fetch-limited 
data described by Ardhui n et alJ (120071) . 

The cumulative breaking term represents the smooth- 
ing of the surface by big breakers with celerity 
C that wipe out smalle r waves of phase speed C 
dBabanin and Youngll2005l) . Due to uncertainties in the 
estimation of this e ffect from observat ions, we use the 
theoretical model of Filip ot et alJ d2008). Briefly, the rel- 
ative velocity of the crests is the norm of the vector differ- 
ence, Ac = |C — C'|, and the dissipation rate of short 
wave is simply the rate of passage of the large breaker 
over short waves, i.e. the integral of AcA(C)dC, where 
A(C)dC is the length of breaking crests per unit surface 
that have velocity components b etween Cx and C x +dC x , 
and between C y and C y + dC y (lPhillipslll985h. Because 
there is no consensus on the form of A dGemmrich et al.l 
2008), we prefer to link A to breaking probabilities. 
Based on Banner et al. (2000, figure 6), and taking their 
saturation parameter e to be of the order of 1.6y/B, the 
breaking probability of dominant waves waves is approx- 
imately P = 28.4 (max{VB - y/B^, 0}j . In this ex- 
pression, a division by 2 was included to account for 
the fact that their breaking probabilities was defined for 
waves detected using a zero-crossing analysis, which un- 
derstimates the number of dominant waves because at any 
given time only one wave is present, and thus low waves 
of the dominant scale are not counted when shorter but 
higher waves are present. 
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Extrapolating this result to higher frequencies, and as- 
suming that the spectral density of crest length per unit 
surface Z(k), in the wavenumber spectral space, is Z(k) = 
l/(27r 2 fc), we define a spectral density of breaking crest 
length, A(k) = Z(k)P(k), giving the source term, 

x max | y/B(f',6' - \/B^, o} ^dd'df 

(B13) 

The tuning coefficient C3 which was expected to be of 
order 1, was here adjusted to 0.4. The resulting model 
results appear to be very accurate for sea states with sig- 
nificant wave heights up to 8 m. Larger wave heights are 
underestimated. Other parameter adjustments can correct 
for this defect, e.g. reducing s u and increasing C3, but 
then the Stokes drift may not be so well reproduced, es- 
pecially for the average conditions discussed here. These 
different possible adjustments and their effects will be 
discussed elsewhere. 



of interest. The only two parameters that are measured 
continuously offshore of the area of interest are the wave 
height H s and mean period /02, recorded at buoy 62163, 
150 km to the west of point A. H s and /02 can be com- 
bined to give the second moment of the wave spectrum 
m 2 = (0.25# s / 02 ) 2 . 

Because there is no reliable wave measurement with 
spectral information in deep water off the French North- 
East Atlantic coast, we also use buoy data and model re- 
sult in a relatively similar wave environment, at the loca- 
tion of buoy 46005, 650 km off Aberdeen (WA), on the 
U.S. Pacific coast. Since this buoy is not directional we 
first examine the third moment of the wave spectrum 

m 3 (/c) = f° fE(f)df. (CI) 
Jo 

If waves were all in the same direction, 7723 would be pro- 
portional to the Stokes drift U ss (fc) of waves with fre- 
quency up to f c , as given by eq. (0). We thus define a 
non-directional Stokes drift 



b. Numerical schemes and model settings 

Spatial advection in the finer model grid is performed 
using the explicit CRD-N scheme (Contour integration 
based Residual Distribution - Narrow stencil scheme 
ICsik et al.l [2002b tha t was applied to the Wave Action 
Equation by Roland (2008) and provided as a module for 
the WWIII model. The scheme is first order in time and 
space, it is conservative and monotone. 

All model grids are forced by 6-hourly wind anal- 
ysis at 0.5 degree resolution, provided by ECMWF. 
The model spectral grid has 24 regularly spaced direc- 
tions, and extends from 0.037 to / max = 0.72 Hz 
with 32 frequencies exponentially spaced. The model 
thus covers the full range of frequencies that con- 
tribute most to the filtered Stokes drift Usf- The 
usual high frequency tail proportional to /~ 5 is only 
imposed for frequencies larger than the diagnostic fre- 
quency fd = Ff m fl,-i, with the mean frequency defined 
by/ m ,o,-i = [/i(/)//d///£(/)d/]- 1 . Here we 
take a factor F = 10, instead of the usual value of 2.5 
dBidlot et al.l 120071) . so that fd is almost always larger 
than the model maximum frequency of 0.72 Hz. Be- 
sides, the time step for integration of the source func- 
tion is adaptatively refined from 150 s for the local model 
down to 10 s if needed, so that virtually no li miter con- 
strains the wave field evolution (Tolman 2002). 

APPENDIX C Model accuracy for relevant 
parameters 

In order to define the errors on the estimations of Us / 
used to determine the quasi-Eulerian velocity U e from 
the radar measurement, it is necessary to examine the 
quality of the wind forcing and model results in the area 



(2ir) 3 m 3 (f c )/g. 
Looking at buoy data we found that 



(C2) 



U S snd{fc) 



5.9 x 10" 



1.25-0.25 



0.5 

To 



1.3 



U 



10 



x min{f7io, 14.5} + 0.027 (H s - 0.4), 

(C3) 

where f c is in Hertz, C/ 10 is in meters per second, and H s 
is in meters. 

Taking directionality into account eq. (0]l yields 
Uss(fc) — 0.85C/ ssn d(/ c ), for typical wave spectra, and 
the relationship dC3b becomes eq. (fT]). For buoy 46005, 
which is a 6 m NOMAD buoy, and f c in the range 0.3 to 
0.5 Hz, this relationship gives a root mean square (r. m. 
s.) error less than 1.0 cm s , corresponding to less than 
15% of the r. m. s. value estimated using eq. ( IC2b . This 
is smaller than the error of estimates using previous wave 
models (24% with the parameterization by Bidlot et al. 
2007), but comparable to the 14.2% error obtained with 
the present model. The same analysis was performed, 
with similar results, for very different sea states recorded 
by NDBC buoys 51001 (North-East of Hawaii), 41002 
(U.S. East coast), 46047 (Tanner Banks, California), and 
42036 (Gulf of Mexico). 

Another source of continous wave measurements is 
provided by altime ter-derived H x , which we correct 
for bias following Oueffeulou (120041) . and fourth spec- 
tral moment 7714. The latter is approximately given by 
dVandemark et alj|2004l) 



TO4 



0.64 5 2 



(C4) 
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FIG. Cl. Variation of the wave spectrum third moment, mg converted 
to a velocity C/ ssn d = (27r) 3 rfi3 (fc)/g, that would equal the surface 
Stokes drift in deep water if all waves propagated in the same direction. 
For each data source a cut-off frequency of f c = fg = 0.36 Hz is 
taken and the data is binned wind speed, at 1 m s — 1 intervals, and sig- 
nificant wave height H s (in colors) at 1 m intervals from 1 to 1 1 m. The 
top panel shows buoy data offshore of Oregon (NDBC buoy 46005), the 
middle pannel shows present model results, and the bottom panel shows 
results from the same model but using the parameterization of Bidlot et 
al. (2007), including a factor F = 2.5. The vertical error bars indicate 
plus and minus half the standard deviation of the data values in each 
(Uio,H a ) class. 



where <tq is the normalized radar cross-section, corrected 
for a 1 .2 dB bias on the C-band altimeter of JASON in or- 



Table Cl. Model accuracy for measured wave parameters in various 
regions of the world ocean. Buoy validation span the entire year 2007, 
except for buoy 62069 for which data covers the time frame 25 Jan- 
uary to 20 August 2008, buoy Iroise covers 13 April to 20 May 2004, 
and JASON 1 data corresponds to January to July 2007 for the global 
validation (JAS-Glo: 393382 data points) and the full year for a box 
3° by 4° centered on 48.5 N and 8 W or 45 N and 128 W. (JAS-Gas 
or JAS-Was: 380 data points). Unless otherwise specified by the num- 
ber in parenthesis, the cut-off frequency is take to be 0.5 Hz, C stands 
for C-band and fg = 0.36 Hz corresponds to our 12 MHz HF radar. 
The normalized bias (NB) is defined as the bias divided by the r.m.s. 
observed value, while the scatter index (SI) is defined as the r.m.s. dif- 
ference between modeled and observed values, after correction for the 
bias, normalized by the r.m.s. observed value, and r is Pearson's cor- 
relation coefficient. Only altimeter data are available at point A but the 
uniform error pattern and the model consistency suggest that errors at A 
should be similar to offshore buoy errors such as found at buoy 62163 
offshore of A, or at the U.S. West coast buoy 46005. Errors at point B, 
not discussed here, are expected to be closer to those at the nearshore 
buoys 62069 and Iroise. 





dataset 


NB(%) 


SI(%) 


r 


2004 


H s 


62163 


6.8 


11.1 


0.977 


fo2 


62163 


10.4 


8.8 


0.907 


H s 


Iroise 


12.8 


17.4 


0.975 


fo2 


Iroise 


-10.0 


11.7 


0.913 


Ussndifs) 


Iroise 


27.2 


26.9 


0.968 


U ss (f B ) 


Iroise 


20.5 


18.5 


0.971 


2007/2008 




JAS-Glo 


-0.6 


11.4 


0.966 


TO 4 (C) 


JAS-Glo 


0.6 


9.1 


0.939 


H s 


62163 


-1.4 


8.8 


0.985 


fo2 


62163 


6.3 


7.3 


0.938 


H s 


62069 


10.1 


14.1 


0.974 




62069 


-7.7 


11.8 


0.886 




62069 


15.8 


24.1 


0.955 


C^ssnd(/s) 


62069 


13.9 


23.0 


0.965 


UssUb) 


62069 


11.1 


21.0 


0.963 


H s 


JAS-Gas 


-2.6 


8.8 


0.983 


m4(C) 


JAS-Gas 


1.0 


6.7 


0.962 


H s 


46005 


4.9 


10.2 


0.975 


f02 


46005 


-2.8 


6.6 


0.931 




46005 


-5.4 


13.5 


0.965 


£4snd( /b) 


46005 


-4.9 


12.6 


0.973 


C/ssnd(0.5) 


46005 


6.2 


12.7 


0.971 


H s 


JAS-Was 


2.4 


7.9 


0.985 


m 4 (C) 


JAS-Was 


1.8 


7.3 


0.953 



der to fit airborne observations (lHauser et alJl2008l) . The 
model estimation of 7714 (0.72 Hz) is extrapolated to C- 
band by the addition of a constant 0.011g 2 /(27r) 4 , con- 
sistent wit h the saturation of the s hort wave slopes ob- 
served by IVandemark et alj (120041) . For this parameter, 
the model is found to be very accurate, especially around 
the region of interest, relatively more so than on the U.S. 
Pacific coast. 

These indirect validations suggest that the third spec- 
tral moment including waves up to the Bragg frequency 
fs = 0.36 Hz, which is proportional to U ssn d, is prob- 
ably estimated with bias between -5 and 10%, and an 
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r.m.s. error less than 20%. The bias on the significant 
wave height appears to increase from offshore (altime- 
ter and buoy 62163 data), to the coast (buoys Iroise and 
62069), and we attribute this effect to the tidal currents, 
not included in the present wave model, and coastal mod- 
ifications of the winds that are not well reproduced at this 
10-20 km scale by the ECMWF model. Because the cho- 
sen area of interest lies offshore of the area where currents 
are strongest (figure 1), we shall assume that, at this site, 
the model bias on U ss (fB) is zero, which appears most 
likely. Extreme biases of ±10% only result in deflections 
of 5 degrees on the diagnosed quasi-Eulerian current Ue- 
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